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ABSTRACT 

We investigate the propagation of electromagnetic fields and potentials in the plasma 
of the early Universe, assuming a Friedmann-Robertson- Walker background with nega- 
tive curvature. Taking over results from classical plasma physics, we show that charged 
particles will acquire an effective mass that has not only the expected thermal com- 
ponent but also a non-thermal component due to the influence of distant matter. 
Although this is a direct effect of the vector potential, we show the theory is nev- 
ertheless gauge invariant. This phenomenon is therefore in the same category as the 
Aharonov-Bohm effect. The non-thermal component becomes increasingly important 
with time, and in some cosmological models can prove to be of decisive importance in 
bringing about the phase transition that generates normal masses. 
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1 INTRODUCTION 



In this paper we study the effective mass of charged parti- 
cles in a plasma, with particular emphasis on the early Uni- 
' verse, before ordinary masses are generated. In Minkowski 
, space, particles of mass m acquire an additional effective 
■ thermal mass, nith, of general order of the temperature of the 
plas ma. This is normally derive d using path-integral meth- 
ods (|Narlikar fc Padmanabhan|[l9 86). We will estimate the 
effective mass from a different point of view, which will al- 
low a straightforward extension to an expanding Friedmann- 
Robertson- Walker (FRW) space. 

Our main conclusion is that if the FRW space has nega- 
tive curvature, the effective mass acquires a new non-thermal 
component, mnt, in addition to the thermal one. This new 
component is due to the cumulative effect of distant matter; 
it is not constant, but increases with time relative to the 
thermal component. In today's standard cosmology, where 
inflation leads to a flat space, mnt will not appear. But there 
are alternative cosmological models in which it could be of 
decisive importance. 

rUnt appears because in curved space the potentials (but 
not the fields) have a tail llDeWitt fc Brehmelll960l : iNarUkaJ 

Il970l : lHaghighipouil[2005l ). as shown in figures [2] and |4] In 
the tail, we have potentials but no fields, so rrint is a direct 
effect of the potentials. This is a quantum mechanical result, 
though we have derived it in the conte xt of classical grav- 
ity. As with the Aharonov-B ohm effect ijAharanov fc Bohml 

ll959l : lTonomura et alllToS^ ). the appearance of mnt at first 



seems to violate gauge invariance, and we have to show why 
this is actually not so. 



2 NOTATION 

We consider a model consisting of a single massless, charged 
scalar field, <E>, interacting with the usual electromagnetic 
fields and potentials. Units are ch osen so that c = h = 1. 
Our sign conventions are those of IWeinberd l| 19721 ). so the 

metric signature in Minkowski space is ( h When 

gravitational effects are introduced, g = — Det {Qfiv)- 

For the early Universe we assume a FRW metric cor- 
responding to negative curvature. (Positive curvature would 
give similar effects, but involves an awkward discrete dis- 
tribution of wave numbers.) We use coordinates t, x, 9, (j>, 
labeled 0, 1, 2, 3, so that 



ds" 



-df^ + R^(t) (dx^ + sinh^ xde^ 
-I- sinh^ X sin^ 6 d(/)^) , 
diag [-l,R^{t),R\t) sinh^ x, 

(t) sinh^ X sin^ 0^ , 
_R^(f)sinh^Xsin6'. 



(1) 

(2) 
(3) 



The expansion parameter, R{t), has the dimension of length, 
and can be thought of as the radius of the Universe at time 
t. 

The conformal time, rj, is related to the ordinary time 
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thy At) = At/R{t). In terms of r), x, d, (j) (s^U dimensionless) 
we have: 



+ sinh^ X sin^ 8 <^4>^^ 



Qfiv = -R^(r?)diag (-1, 1, sinh^ XjSinh^ Xsin^ 
^/g = i?'*(77) sinh^ X sin ^• 

Maxwell's equations in free space are: 



al3\ 



= 0, 
= 0. 



(4) 
(5) 
(6) 



(7) 

(8) 



Potentials and a gauge condition will be introduced 



later. 



3 EFFECTIVE THERMAL MASS IN A 
PLASMA IN MINKOWSKI SPACE 

We start from the simple, and at first sight pointless, ob- 
servation that the kinetic-energy term of the Klein-Gordon 
equation for a charged scalar field in Minkowski space, 
^*(p — when expanded, gives the term $*g^A^$. 

This has the same sign as the mass term "l>*m^"I>, suggest- 
ing that in some circumstances, such as in a plasma, {q^A?) 
might play the role of the square of an effective thermal 
mass. (A^) is the sum of contributions from particles in 
the immedi ate neighb orhood, roughly, those within the De- 
bye sphere l|Sturrocklil994l ') . and should therefore be written 
((^. Ai)^). When we sum over all contributions, the Ai po- 
tentials will be randomly orien ted, and will add according 
to the theory of random flights jHughedfiggSf ). Ai) wfll 
then be zero, but {(^ Ai)^) will not. 

In a proper quantum mechanical treatment we have to 
represent the source particles by means of a density matrix. 
This is not important for a plasma in Minkowski space, but 
will be essential when we come to curved space, so we intro- 
duce the idea here. The density matrix provides a way of in- 
cluding all possible configurations of source particles, which 
are not, of course, individually observed. Instead of Ai we 
should write A^i, where ^ specifies a particular configura- 
tion. (A'^) is then a double sum, over the particles within a 
configuration, and over the various possible configurations; 
up to a normalization factor, (^^(^^ A^i)^). Dimensional 
considerations indicate that in the conditions of the early 
Universe this should have a value of order T'^, where T is 
the temperature of the plasma. 

We are free to make an overall gauge transformation 
with an arbitrary function i/;(a;^). The various Ai become 
Ai -I- V?/', and the source particles acquire corresponding 
phases. But we should note that this freedom is not the 
same as freedom to choose any gauge we wish to describe 
the propagation of the potentials. For example, we might de- 
cide to switch from Lorenz gauge (jjacksonl 19991 ) to Coulomb 
gauge. This would introduce a separate gauge function for 
each A^i, and there would be no single gauge function that 
could be used to generate a compensating phase for any 
chosen interacting particle. If we are concerned with the 
propagation of potentials we have to choose a gauge from 
the beginning that is suitable for that purpose, as discussed 



in section [T] This choice is not altered by an overall gauge 
transformation. 

We can argue, of course, that in Minkowski space we are 
not really dealing with a direct effect of the vector potential, 
because electric and magnetic fields are also present, and 
they are the actual physical agents. This is true, but, as we 
will see in section (5] it is no longer true in a curved space, 
where we have potentials without fields. 

We will need to know the potential due to a single 
particle in the plasma at distances larger than the De- 
bye length. If the particle is stationary, the potential de- 
clines exponentially. Sur prisingly, however, starting in the 
1950 ' s several authors jNeufeld fc Ritchiel 1 19551 : iTappertI 
1 19671 : iMontgomerv et al.ll 19681 ) discovered that a test charge 
moving uniformly is not exponentially screened, but gener- 
ates the field of a quadrupole in a collisionless plasma. The 
restriction to a c oUisionless plasm a was lifted in later papers 
(|Yu et al.lll973 : ISchroeder|[l975l ). where it was shown that 
the far field of a moving test charge is that of a dipole with 
strength of order qrtV , where q is the charge, Tt the colli- 
sion time as measured by t, and V the velocity. This dipole 
form ensures that in Minkowski space the contribution of 
particles to the effective mass falls off rapidly with distance 
beyond the Debye radius. We shall see, however, that this is 
no longer the case in a curved space. 



4 SOURCE PARTICLES IN THE PLASMA 

To study the propagation of fields and potentials in a curved 
FRW space, we first need a model of the motion of the 
source particles that produce them. We are concerned with 
the early Universe, before any symmetry-breaking transi- 
tions that may generate particle masses as we know them. 
The properties of the plasma are therefore determined by a 
single number, the temperature. The photons will have the 
familiar Planck distribution. For simplicity, we assume the 
charged particles are associated with a single massless scalar 
field. This field will have an effective thermal mass, mth, due 
to its interaction with the photons, so that mth ~ T. The 
charged particles will have a thermal distribution appropri- 
ate to this mass and temperature. The plasma will there- 
fore be weakly relativistic; for simplicity, however, we will 
carry over some results from non-relativistic plasmas. Other 
plasma parameters are determined by T; for example, the 
plasma frequency, ujp, will be of order T. 

A charge that is part of the plasma, rather than being 
a test charge constrained to move in a certain way, will it- 
self experience collisions. We can picture such a particle as 
describing Brownian motion, with its screening cloud some- 
times closer, sometimes farther, but never fully established. 
We will model the most important aspect of a single step 
of this process by setting up a local set of Cartesian co- 
ordinates and placing a stationary screening charge — q at 
the origin. A charge q is initially at rest on the z axis at 
z — —Vrt, then begins to move with uniform velocity V 
along the z axis, starting ai t — —Tt and ending at t — Tt, 
when it again comes to rest. Note that this choice of origin 
for t differs from that of the cosmic time used in section [2] 
rj, defined below, will differ in a similar way. We reconcile 
these differences later, in section [TO] 

We are now going to take the result for a uniformly 
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moving test charge and use it to calculate the field around 
our model charge that is subject to Brownian motion. We 
justify this step as follows. Imagine that the test charge, 
instead of moving uniformly from an indefinite time in the 
past, is actually stationary until time t = 0, and then starts 
moving uniformly. There will be transient fields, but after 
a few collision times the final field of a moving dipole will 
be established. During the transient period, the test charge 
will be moving away from its screening charge, which only 
gradually picks up speed and trails along behind. This is 
similar to the motion of the charge in our model. Since the 
transient fields form a bridge between the initial and final 
fields they must have a similar form to the field of the test 
charge, and extend comparably far. 

During the period that the charge is moving (the pulse), 
the dipole moment will be D{t) = qVt. Before and after the 
pulse the dipole moment remains at a constant value, but 
this is of little interest because a constant dipole moment 
generates no magnetic field or vector potential. It is well 
known that Maxwell's equations separate in conformal time, 
so we will write the dipole moment in terms of rj. Define 
r = Tt/Rs{t) (dimensionless), so that the pulse extends from 
rj = — T to 7; = r. The subscript 's' on Rs indicates the source 
time. Rs{t) can be treated as constant during the pulse, and 
the dipole moment can be expressed as a Fourier integral: 



^^^^ - h 



D{r, 



D{n) exp(— in?;) An, 
2iqVRs sin(nr) 



(9) 
(10) 



We estimate the value of rt in section fTUl The z axis of 
the local system used in this section will be parallel to the 
polar axis of the polar coordinates used in the bulk of the 
paper. 



5 DIPOLE FIELDS 

We will need only the three fields _Fio, -F20 and F12. Since 
Maxwell's equations separate in conformal time the elemen- 
tary solutions can be written 



Fio = /io(x,^',0)exp(-inr7), 
F20 = /2o(x,^',0)exp(-inr7), 
Fx2 = /i2(x,^',0)exp(-in?7). 



(11) 
(12) 
(13) 



For dipole fields we try the following forms: 
/lo = h{x)Pi{0), (14) 



/20 = f2ix)Ni{e), 
fi2 = Mx)Ni{e), 



(15) 
(16) 



with angular functions Pi (6) — cos 9 and A''i (6) = 
dPi(6i)/d6l = - sine. 

Maxwell's equations now give (with a prime meaning 
d/dx): 



insinh x/i ^ 2/3 = 0, 
- = 0, 
h-f2+^nh = 0. 



(17) 
(18) 
(19) 



5.1 Dipole Magnetic Field 

From (|17p . (|18|) and (|19p we obtain the equation for /a alone 
2 



dVa , / 2 



/3 = 0. 



(20) 



dx^ ' V " sinh^x. 

This is the analog of equation 16 of iMashhoonI l)l973h . 
which he derived for a space of positive curvature. Define 
Uix) = faix)/ sinhx- Then satisfies 



(21) 



This is a special case, for / = 1 or Z = —2, of th e equation for 
hyperbolic spherical functions, with solution l|Bucher et al.l 
1 19951 : iBander fc Itzvksonlll966l ) : 

/4(/,x) =iVK-l)'+^sinh'x^(^^|J^cos(nx) (22) 

where A*'; is a normalization factor that is irrelevant here. 
Choosing I = 1 (—2) we get the solutions of (|2H) that are 
regular (singular) at x = 0. With / + 1 < interpreted as 
integration: 

/4,rog = /4(l,x) = 3-J2 — ["-sinhx cos (nx) 



sinh X 
— cosh X sin(nx)] 



/4,sing = /4( — 2, x) — 



[n sinhxsin(nx) 



sinh^ X ' 
+ cosh X cos(nx)] 



(23) 



(24) 



/a is constructed from that combination of /4,rog and 
/4,sing that is proportional to exp(mx), because when com- 
bined with exp(— in?;) this represents outgoing waves: 



. Ci{n) exp(inx) , , .... 
fa = — , — — (cosh y — m smh x) , 
smh X 



(25) 



where Ci (n) is a normalizing factor to be determined. 

It is convenient at this point to express /a in terms of 
u — tanh(x/2): 

Ci(n)exp(inx) 



2u 



(1 — 2in-u + 



(26) 



5.2 Dipole Electric Fields 



From (|17p and (I25p we derive the equation for the radial 
electric dipole field: 

^ -iC,in)e.pJ^nx)il-ur _ ^.^^ ^ ^2) ^^^^ 

Similarly, from (I18|l and (|26p we derive the equation for 
the transverse electric field: 

iC'i(w)exp(inx) 
72 = : X 



[(1 - u^f - 2inu{l + u^) - Anu\ 



(28) 



5.3 Normalizing the Fields 

From ([9| and pO|) we can get the near-field expression for 
the radial component of E, and by comparison with (|27|l 
arrive at the form of the normalizing factor Ci(n). On the 
polar axis, at small distances, the radial E field at frequency 
n is 
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Figure 1. Typical contour for Fourier synthesis of magnetic and 
electric fields. 



^ ^ d_ f Djn) exp(~m??) \ ^ 2D{n) expj-mr]) 

^ dr \ r'^ J 

where r = -RsX- 

In the same hmit (small x)j l|27p gives /i(x) = 
— 2iCi (n)/(nx^), and so, on the axis, 



^'io(r?,x) = -2iCi(n)i?f exp(-irz77)/(nr'). 



(30) 



We now transform to a local Minkowski frame with co- 
ordinates t, r: 



= — 2iCi(n)-Rs exp(^in7;)/(nr^). 



(31) 



-Fio is just the conventional radial electric field as given 
in ((m, so using (fTO)) : 



Ci(n) = 



~2qV sin(nr) 



(32) 



6 PROPAGATION OF THE MAGNETIC 
FIELD 



(|26|) . (|32|) . (fT6|) and ([T3|) give the Fourier transform of the 
magnetic field: 

--(7Vsin(nr)iVi(6') exp[in(x - ??)] 
ri2 = X 



(l — 2inii + 



(33) 



We transform back into x, V space by dividing by 27r 
and integrating over n along the real axis. For x > '7 + ''" both 
exponentials in sin(nr) allow us to close in the UHP, using 
the contour of figure [T] sin(nr)/n is regular at n = 0, so we 
get zero, as required by causality. Similarly, when X ^ V~ 
we can close in the LHP, using a contour that is the inverse 
of figure m and we again get zero. 

For rj — T < X < ??-l-rwe must write sin(nr) = 
(exp(inr) — exp(— inr))/2i, divide the integrand into two 
pieces, and close the first integral in the UHP and the sec- 
ond in the LHP. Combining these two integrals we get 



Fl2 



qV smB{\ + \ 
2^ 



(34) 



The propagation of -F12 is shown in figure [2] An unex- 
pected feature of the pulse is that the amplitude does not 
tend to zero for large x, but to a constant asymptotic value. 
The conventional H field, does, of course, tend to zero, in 
fact exponentially. We can see this at a particular instant 
by setting up local Cartesian coordinates at a particular 
observation point, P, at rest in the FRW space. Choose a 
point in the equatorial plane, Q = 7r/2, with azimuth (ji = 0. 



I I 



I I 



I I 



X-O 1 2 3 4 5 6 

Figure 2. Propagation of the magnetic field, -F12: snapshots of 
the pulse for various values of conformal time, r\. For each 77, the 
pulse is confined to the range r] — T < x < '? + ''"; within that range 
we plot (1 + v?)/u. To render the pulse visible on the graph, we 
have arbitrarily set t = 0.5, far larger than it is in practice. 



Choose dx parallel to dx, dy parallel to AO, and Az paral- 
lel to A(j). Then da; = -Rpdx, and Ay — i?psinhxd6, where 
the subscript 'p' refers to the particular instant chosen. The 
conventional field is F12, given by 



F12 



dx dy 



(35) 



= Fi2/(i?^sinhx), 
showing the exponential decline for large x- 

7 DIPOLE POTENTIALS 

The potentials are defined by the following equations: 

Ao = ho{x)Pi{e)exv{-iniq), (36) 
Ai = /ii(x)Pi(e)exp(-inJ7), (37) 
A2 = h2ix)Niie)exp{~mrj), (38) 
dA^ _ dA^ 
dxt^ dx'^ 
From (1391) we derive: 



(39) 



^2 — h\ 
ho + vnhi 
ho + mh2 



fi, 



(40) 
(41) 
(42) 



(l40l) . (|4T]) and (|42]) are not independent. If we differ- 
entiate (|42p . subtract it from (|4ip and use p9p . we obtain 
(1401) . These three equations therefore do not suffice to de- 
fine the potentials uniquely; this is to be expected, because 
Af_i is only defined up to a scalar gauge function, 1/;, so that 
Afi — Afi -\- dij^/dx^ defines the same fields as j4,t. 

To fix the potentials uniquely we need a gauge con- 
dition. Although, as we saw in section [31 there is a real 
sense in which our final result is gauge invariant, we still 
have to specify a gauge condition at this point to determine 
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how the potentials in our model propagate over cosmolog- 
ical distances. T he natural choi ce is the Lorenz condition 
(LC), as given in IjacksonI ()l999l ). This is covariant in form 
and respects causality; it also has the advantage that the 
current density generating the vector potential is the com- 
plete current density, not (as for the Coulomb gauge) just 
the transverse part. The curved-space generalization of the 
LC is A'^,, — 0, which in our metric becomes 



0. 



(43) 



This has the disadvantage, in a FRW space, that fac- 
tors of R{ri) remain. This means that the equations describ- 
ing the propagation of the potentials will involve the time 
derivative of R, and so will be coupled to the gravitational 
field equations. It is certainly possible that Nature works 
this way. But to simplify the mathematics we propose here 
to use the following modified Lorenz condition (MLC): 

^ {V9 9'"'A./R')^0, (44) 



dxt^ 
which leads to 



insinh^ xho H fsinh^ X^^i] " 2/i2 — 0. 

dx 



(45) 



(1441) may seem unphysical, because it depends explic- 
itly on R{rj), and so seems to be restricted to FRW spaces, 
rather than being universally applicable. But it can be given 
a physical meaning in the context of conformally invariant 
theories. The connection is outlined in appendix [B] 

The potentials are normalized by comparing them with 
the expressions for the fields. Together with the MLC, this 
determines the potentials uniquely. 



T.l Scalar potential 

By combining HT]), ^8^, (gl]), ((42]) and (gS]) we can derive 
an equation for ho alone: 



-— -+2cothx-j— + n - . ,2 
dx dx V smh-'x 



ho = 0. 



(46) 



Defining a = \/n? — 1, this becomes the equation for 
hyperbolic spherical functions, (|2H) . As before, we choose 
the solution that represents outgoing waves. With u — 
tanh(x/2): 

Ao = Pi{B)ho{x) exp{—mri) 

= C2 (n)Pi {6) exp[i(ax - n-q)] X 
(1 - u'^) (1 - 2iau + u^) 

-\-2 (47) 

where C2{n) is a normalization factor that depends on the 
driving function. C2(n) is, of course, proportional to the 
normalization factor Ci(n) defined in (|32p . For small x, (147p . 
1(42)) and ^ give: 



(48) 



C2(n) = 



iCi(n) 
4n 



7.2 Vector Potential 

The interesting component of the vector potential is the 
transverse one, Ae, or in our notation A2. (|42p gives 



Figure 3. Integration contour for A2. The contour encloses the 
points n = —1 and n = +1; it is described in a clockwise sense. 



(49) 



This can be used to get /12, and from that A2, since /2 and 
ho are known from (|28|) and (f47|) . 

Equation (|49|) shows that A2 is the sum of two parts, 
one involving exp(inx), the other involving exp(iax)- 



A2 ^ -a^MWsin(n^)e— (^^^)^ 
2n^ij^ 



(50) 



N = e'''''[{l-u''f ~2mu{l+u^) - ^r?u^] , (51) 
A = e'°"'{l-u^)[-l + 2iau-u^] . (52) 



8 PROPAGATION OF THE VECTOR 
POTENTIAL 

The Fourier synthesis of A2 is carried out, as for the mag- 
netic field, by dividing by 2-k and then integrating over n, 
along a contour chosen to respect causality. There is no vec- 
tor potential before the pulse begins, so the correct contour 
is a line parallel to the real axis and slightly above it. The 
integral must be taken along the whole path A-D in figure[TJ 
including the semicircle B-C. This ensures that for x > 77 + 1" 
the contour can be closed in the UHP and the integral will 
be zero. There are two regions of interest, ry — r<x<'7 + ''" 
(the "main pulse"), and X < V ~ (the "tail"). 

We get a non-zero result only for that part of the inte- 
gral that involves a contour that is closed in the LHP. For 
the part of equation H50|) that is proportional to exp(inx) 
we can shrink this contour to a small circle about the point 
n — 0, and we just have to find the residue there. For the 
part that is proportional to exp(iax) have to remember 
the branch points at n = ±1, so our contour can only be 
shrunk to the form shown in figure O 

A2 is the sum of contributions from small circles (or 
semicircles) around n = (the pole terms) and the inte- 
grals along the cut from n = — lton=l. The pole terms 
have simple analytic expressions, but the integrals must be 
evaluated numerically for each {rj, x) pair. The integrations 
are straightforward, and the propagation of A2 is shown in 
figure m 

An important difference between figures [2] and [4] is that 
in the latter the pulse has a non-zero tail for x < rj — r. This 
feature of the propaga tion of potentials in a cur v ed space 
has been noted before llDeWitt fc Brehmd Il960l): iNarUkaj 



HaghighipouJ (|2005l )). For 



(l973); for more recent work seel 
dipole propagation, as here, the tail rises linearly from x = 0, 
and approaches a constant value for large X- This asymp- 
totic value is proportional to r. In ordinary Minkowski space, 
which corresponds to the limit x ^ 0, A has no tail. 

We note that the following simple function, with a = 
0.735, gives an adequate fit for A2 for large 77: 



Ai 



(2gVr sin 9)^ tanh^ (ax) . 



(53) 
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CALCULATION 

We will show here that (A'^A^j,) is composed of a thermal 
part, {j4''ylfi)th, which will be of order T^, and a non-thermal 
part, {A'^A^} nt- We are concerned with the ratio 





X-O 1 2 3 4 5 6 

Figure 4. Propagation of the transverse component of the vector 
potential, A2: snapshots of the pulse and the tail for a series of 
values of rj. As with the magnetic field, we set r = 0.5 to make 
the pulse visible on the graph. 



9 EFFECTIVE THERMAL MASS IN A 
PLASMA IN FRW SPACE 

In section [3] we discussed the thermal mass in Minkowski 
space; we now have to investigate any changes there may be 
when we go to a FRW space. The most important one is the 
appearance of the tail following the pulses of A. Since the 
fields themselves do not have tails, the vector potential in 
this region is a pure potential, without accompanying fields. 
It can therefore be represented as the gradient of a scalar 
function, dS/dx'^. 

In the asymptotic region, where both rj ^ 1 and x ^ 1, 
S has a simple form. Both ho and hi tend to zero exponen- 
tially as X ~* 00 ■ ''■2 is given by the contour integral of the 
previous section, and in the asymptotic region the integral 
along the cut is negligible compared to the pole term. So we 
get A2 ~ 2qVTNi{e), and 5 ^ 2qVTPi(e). 

If the potential from source particle i is denoted by A2i, 
then the total from all the particles in any one configuration 
is iS = ^ . iSi . We might think that we could use 5 as a gauge 
function and remove the effect of A2 in the tail entirely. 
But this is not so once we introduce a density matrix and 
sum over all configurations of the distant source particles. 
Each configuration now produces its own gauge function, 
Sfj, = 5pi. These functions are all distinct, and there is 
no single gauge function that can be used to eliminate the 
effect of A2 in the tail. 

This is not to say gauge invariance is violated; as dis- 
cussed in section [31 we can introduce an arbitrary overall 
gauge function %l){x'^). But this is quite different from the 
separate functions S^i, each of which is determined by the 
position and orientation of the corresponding source dipole. 



P={A^A^)nt/{A^A^)a,. 



(54) 



and its variation with time. 

Moving particles in the distant plasma will generate 
pulses of Af_i, each of which consists of a "main pulse" and a 
"tail" . The main pulse contains electric and magnetic fields, 
which are distinguished from the thermal background by 
their spectrum. The associated vector potential will con- 
tribute to {A^A^)nt- We will not investigate this in detail, 
because, as we will see later, it is the potentials in the tail, 
not the main pulse, that are responsible for the increase of 
p with time. A^ in the tail is a pure gauge potential that 
generates no fields, but it, too, will contribute to {A^A^)nt- 

To study the development of {A'^A^)nt we need a model 
of how the Universe came into being. We will s uppose it orig- 
inate d as a bubble in a metastable vacuum ( [Bucher et al.l 
Il995h . and for simplicity assume that all fields at tempera- 
ture Tinax appeared without appreciable delay at the surface 
of the bubble. The evolution of {A'^Afj,)nt can be visual- 
ized as follows. As each pulse passes the observation point 
it leaves a memory in the form of the tail. These tails will 
not ca ncel but will a dd according to the theory of random 
flights (|Hughej|l995l '). The scalar potential, of course, will 
remain close to zero; it is only the vector potential that accu- 
mulates these additions. Consequently, from now on we will 
write {A'^An)th = (A^)th, and similarly for the non-thermal 
part. 

In (|54p . both the numerator and the denominator will 
decrease as the temperature falls, but the numerator falls 
more slowly, so the ratio will gradually increase from zero. 
In some cosmological models, as we will see in section II6I 
we would like to know when p becomes of order unity, so we 
will estimate this here. 

It is convenient at this point to change coordinates so 
the observation point is at 77 = rjp, x = 0, and a general 
source point is at rj, x- Our formulae involve differences in 
rj and x, which are unchanged by this shift of origin. We 
will consider ry and t as starting from zero at the time of 
minimum radius, -Rmin, and maximum temperature, Tmax- 
We can picture the buildup of {A^)nt with the help of fig- 
ure [5] where the coordinates are rj (upwards) and x- The 
row of boxes at the bottom of the diagram is at 77 = 0; this 
represents the surface of the bubble. Each box has thick- 
ness dx and duration r. Two observation points are shown, 
Pi and P2, at times 771 and 772. The lines Pi-B and P2-C 
represent the past light cones. All plasma particles within 
these past light cones contribute to {A^)nt. The contribu- 
tions add incoherently over times longer than the collision 
time, r, so we can imagine the whole diagram divided into 
time slices of duration r, just like the one shown for 77 = 0. A 
box at source time 77s represents a spherical shell of volume 
4tvRs sinh^ xdx- 

The number density of particles in this shell, for a ther- 
mal distribution, is 
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= 0.24Zf, 



(55) 
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Figure 5. Buildup of (A^)nf Each box represents a splierical 
sliell of radius Xi tliickness dx, duration r. For clarity, only one 
time slice is shown; the plasma actually fills the whole volume 
from J7 = to the observation time. 



and the number of particles in the shell is then 

diVs = 47r(0.24ri')i?,^sinh^Xdx 

= 47r(0.24)y4^/*sinh^XdX- (56) 

Here we have introduced the quantity Am = R^T*^, which 
(at our current level of calculation) will be constant during 
the expansion, and therefore does not need another sub- 
script, 's' or 'p'. 

Multiplying (|56p and (|53p we get the contribution to A2 
from the particles in the shell. We convert this to g'^^ A^^Av 
by multiplying by g^^: 

4.(0.24) A^/^ sinh^ X dx {2qrf {V'){sin' 9) ^^^^ . (57) 

Hp smh X 

The average of sin^ 6 over a sphere gives 2/3. For the 
average of we reason as follows: if the effective of the 
particles were purely thermal, we could write — 1/2. But 
the effective mass actually increases with time as the ratio 
p from (|54[) increases from zero towards one. So a better 
estimate of is 1/(2 +p). 

= at the surface of the bubble, and = % at the ob- 
servation point. For some time slice at time rj between these 
two limits, we get the contribution of all spherical shells 
to (A'^)nt by integrating from x = to the light cone at 
X = Vp-V- 



Ma24) 3/4 
3 



1 



2+p 



{2qry 



"P"", tanh2(ax) 

dx — m — ■ 

Up 



(58) 



We now have to sum over all time slices of duration r. 
The sum can be converted to an integral by multiplying by 
drj/r: 



tanh"(ax) 

dx p2 
-ftp 



(59) 



In section 13] we defined r = Tt/Ra(t). For Tt we will sim - 
ply use equation 3.24 from iMontgomerv fc TidmanI l|l964h . 

The factor in the numerator will be of order {2 + p)~^^^; 



the quantity in the bracket in the denominator is of order 
unity and will be omitted. We then have, in our notation, 
Tt = m'^{2+p)~^^^ /{S-KUsq'^ In A). For we will use the ef- 
fective mass at time t. The purely thermal value is q^T^ , but 
we should multiply this by (l-|-p) to include the non-thermal 
part also. We have Va = 0.24Ts^, so 

l+p 



8n lnA(0.24)Tsg2i?,(2+p)3/2 

1 +p 

8tt In A(0.24)g2 Ajj^''(2 + p)3/2 ' 



(60) 



We can now get an expression for p by using (|60|) in 
((59)) and dividing by Tp-. 



l+p 



V?2<*(2+p)5/2 



tanh2(ax) 
dx- 



RpTp 



drj- 



l+p 



31nA /„ -'(2+^)5/2 

* l+p 
31nA ''(2+p)5/2 ^ 

{Vp - tanh[a(?7p - rj)]/a} . 



dxtanh (ax) 



(61) 



This integral equation is surely not correct in detail, 
but it does provide a general picture of the buildup of p 
from p = at r^p = 0. We seek »7cw, the value of rjp for which 
p — 1. Use a = 0.735, as before, and for definiteness set 
InA = 10; numerical solution then gives r^cw ~ 11. 

A remarkable feature of (|6Hl is that not only has Am 
disappeared, but also the coupling constant, q. The long- 
range character of the Coulomb interaction, however, is still 
apparent in In A. An equation of this sort might reasonably 
have been expected to yield a value of rjcv that was either 
extremely small, so the CW transition takes place almost 
immediately after the formation of the bubble, or extremely 
large, so the transition never takes place at all. Instead we 
get a value of 77cw that is within one or two orders of mag- 
nitude of unity. 



11 INTRODUCTION OF CONDUCTIVITY 

Up to this point we have been mainly concerned with po- 
tentials in the tail. But we have now to look more closely at 
the electric and magnetic fields in the main pulse for dipole 
propagation. These combine to give a Poynting vector di- 
rected outwards. In fiat space, this vector falls off like 1/x^, 
which tends to zero even when integrated over the whole 
sphere. When we go to curved space, however, the Poynt- 
ing vector tends to {qV sm9/ sinhx)'^- The surface area of a 
sphere is 47r_R^ sinh'^ x> so the integral of the Poynting vector 
tends to a constant non-zero value. One effect of curvature 
is to continuously increase the energy of the plasma. 

It seems unlikely this extra energy will remain in the 
form of pulses like those of figure [S] We know such pulses 
propagate unchanged in a fiat space, but when they have 
traveled cosmological distances, so the curvature becomes 
noticeable, we should expect a gradual thermalization. We 
will model this thermalization by including a small, constant 
conductivity, a, in our equations. Note that a is not directly 
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related to the normal conductivity of the plasma, which (in 
our cosmological units) would be very high, a is just a device 
to represent the slow thermalization, and will have a value of 
order unity. For technical reasons, which we explain below, 
we choose a = l/{2n). 

There are two reasons why we have to investigate the 
effect of cr: 

(i) The thermalization of the pulses raises the tempera- 
ture of the plasma, and so tends to reduce the value of the 
ratio p. 

(ii) The slow decline in the height of the pulses will prob- 
ably reduce the value of A in the tail, and this will also 
reduce the value of p. 

The gradual transfer of energy to the plasma implies 
that the product R(rj)T{ri) will not remain constant, as in a 
simple expansion, but will slowly increase. In this respect cr 
simulates inflation, but on a much longer time scale. 

The equations involved here are analogous to those used 
previously, but are considerably more complicated, so we 
just give the main results. 



12 INCLUSION OF CONDUCTIVITY: 
FIELDS 

The wave number, k, the permittivity, e(n), and the conduc- 
tivity, (J, satisfy the dispersion relation 

— n^e{n) = n(n -I- 4nia). (62) 

By analogy with Maxwell's equations in flat space, the 
equations p7|l . p8|) and (|19p for our dipole fields become: 



inesinh x/i ^ 2/3 = 0, 
- /a = 0, 
/i-/2'+in/3 = 0. 



(63) 
(64) 
(65) 



e 



( 



B 



Figure 6. Contour for computing the magnetic field, -F12, when 
conductivity is included. The points on the negative imaginary 
axis labeled A and B are at n = — STricr/S and n = —Airia, re- 
spectively. The contour is traversed in a clockwise sense. 
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12.1 Magnetic Field 

From these we obtain an equation for /s alone: 
2 



sinh X 



/3 =0. 



(66) 



Writing k^ for n^e this equation takes a familiar form; 
/a can be derived from the empty-space formula simply by 
writing k for n everywhere, except, of course, for the nor- 
malization function, which we denote now by C3(n, a): 

/3 = g3(n,a)exp(ifcx) „ ^.^^ ^ ^2. ^g^^ 

C3{n,a) is most easily determined in flat space; this is 
sufficient since we only have to consider small distances. The 
calculation is done in appendix [X] We find: 



Cain, a) = 



-6qV{n + Ania) sin(nr) 
n(3n + STria) 



(68) 



12.2 Propagation of the Magnetic Field 

The Fourier transform of the magnetic field becomes 

—3qV{n + 47ria) A'"i (9) sin(nr) exp(ifcx — inrj) 



{3n + 8'ivia)nu 



Figure 7. Propagation of the magnetic field, -F12, in the presence 
of conductivity, cr = 1/2tt. 



(1 - 2ifcu + u^) 



(69) 



We can integrate around the pole at n = in the same 
way as before, except that we have to respect the branch 
points of fe, at n = and n = —Ama. We have also to take 
account of the pole at 71 = — STria/S. A suitable contour is 
shown in figure [6] 

The integrals are straightforward, and the resulting 
propagation of F\2 is shown in figure [T] 



13 INCLUSION OF CONDUCTIVITY: 
POTENTIALS 

Her e we follow the prescription of Landau and Lifshitz 
(Landau et al.1 ,1954'): use the same equations relating po- 
tentials and fields as in empty space, i.e. (|40[) . (I41|l and (|42|) . 
but modify the Lorenz condition by including the permit- 
tivity in the ho term: 



inesinh^ yhtj -\- -rr— [sinh^ x^^i] 
dx 



2/12 = 0. 



(70) 
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13.1 Scalar Potential 

From these equations, as before, we can derive an equation 
for ho alone: 



J- + 2cotlix^j7- + 



sinli X 



ho = 0. 



(71) 



Tliis is tlie same equation as we obtained for empty space, 
except that, just as for the magnetic field, in place of 
we must write = n^e. (|47p applies as before, provided 
we write a = \/k^ — 1, and use a normalization function 
C'4{n,a) in place of C2{n): 



C4(n, a) = 



iCsjn, a) 
4{n + 4-Kia) ' 



(72) 



13.2 Vector Potential 

Defining Up = — Svria/S, we find: 
-gViVi(6l)sin(nr)e-"" 



2(n — np)Ti?u'^ 
Aa = e'"''(l-it^) (-l + 2iait-it") . 



(73) 

(74) 
(75) 



We note that going from (f50|) to (ff3|) takes three simple 



steps: 

(i) Rename M and A; call them /Co- and Aa, respectively. 

(ii) In the final parenthesis of (|50p . substitute k for n in 
both fCa and Aa\ this includes redefining a to be ^/k^ — 1 
rather than \/n^ ~1. 

(iii) In the prefactor of (|50p . change v? in the denomina- 
tor to n^(n — Tip). 

(|73|) becomes (|50|) in the limit cr ^ 0, as it should. 



13.3 Propagation of the Vector Potential 

For the Fourier synthesis we need, in general, to use a more 
complicated contour than the one shown in figureO because 
we have to take account of the branch points of both k and 
a. With our choice of cr = l/27r, however, the branch points 
of a coalesce into a single point at n = — i, and the contour 
of figure O is adequate. 

The result is shown in figure [S] 



13.4 Asymptotic Form of the Vector Potential 

We can show easily that in the limit -q oo, A2 satisfies a 
diffusion equation. The main terms come from the integral 
down the cut, and in this limit the integrand is concentrated 
close to n = 0. Note that this integral involves only KLa, the 
"field term"; the Aa ("scalar potential") term, being even 
in n, does not contribute. From this integral we obtain a 
simple analytic expression for the asymptotic form: 

3g1^r sin 9 



Ficy,X,v) 



16n 



Fi^,X,v) 
erf I X 



[(1 + M^)?7 + 47r(Titx] exp 



(76) 



(77) 





X-O 1 2 3 4 5 6 

Figure 8. Propagation of the vector potential, A2, in the presence 
of conductivity, a = 1/2it. 




256 



f - 



Figure 9. y'r;A2 plotted against g = x/%/Vy to demonstrate the 
asymptotic form. In this graph we use the asymptotic form l)76p 
for the tail. 



When X is also large, of order y^, the term involving 
the error function is negligible compared to the others. 

In figure |9] we plot ^JrjA^ against ^ = xl \f^-. with A2 
computed using the asymptotic formula. We see ^/rjA^ tend- 
ing to a constant form; with our choice of abscissa, this is 
a Gaussian for large x and 77. For 77 ^ 4, this graph is in- 
distinguishable from the one using integration around the 
contour. 



14 THE ENERGY-TRANSFER EQUATION 

We will now investigate the effect of a non-zero co nductivity 
on the temperature of the plasma. We start from I Weinberg! 
(1972), equation 4.7.9: 



(78) 
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This expression is zero for any value of v, but the most 
important one for us is = 0. We treat the plasma as a 
perfect fluid at rest in FRW coordinates r), x, &, (p, so that the 
4-velocity ((7°, U') has W = 0, i = 1, 2, 3. Since g^.^U^'V = 
-1 and 300 = R^(r)), U° = 1/R. 

Multiply (|78|l by ,/g d^x U^, sum over v, and set the 
resulting expression equal to zero; we get 

^d\UoT''"^^0. (79) 

The left-hand side of this equation is a coordinate scalar 
with y/gd^a; an element of proper volume. Expanding (|79|) : 



V^d xUo 



^gd^A^gT )+— — (^T 



+ 



V°,T'' + rlT""] = 0. 



(80) 



From this we can derive an integrodifferential equation 
for T^". The resulting equation is most simply written in 
terms of j4m(??) = {ri)T'^ {rj) (eV)* • m*. A similar quantity 
was introduced in section 1101 but there it could be treated 
simply as a constant. Here we calculate its variation with 
time: 



Am{v) 



a(77^_^^ 

2+p{r!) 

exp [— 27r(j(r; — 77')] 
P4 = 3.97 X 10"'^ eV • m 



= 0, 



(81) 
(82) 



Note that because of the appearance of the function 
p(?7), this equation by itself is not sufficient to calculate Am- 
We obtain an independent equation connecting these two 
functions in the next section. 



15 MASS GENERATION: 
CALCULATION 



SECOND 



We will now obtain an equation for the ratio p{ri) analogous 
to (|6ip . but based on the asymptotic form (|76p for the tail of 
the vector potential. In our previous calculation, where the 
tail of A2 reached a constant asymptotic value, we could be 
sure that p would rise monotonically; the only question was 
whether there would be factors of Am that would make the 
increase far too fast or far too slow. (This is just what hap- 
pens, for example, if we use quadrupole rather than dipole 
potentials.) In this second calculation, we can anticipate that 
there will not be any unwanted factors of Am, but only a 
detailed calculation will show whether p continues to rise 
with 77. 

The analog of (|56p is 

diV, = 4^(0.24) A^/^ (77s) sinh^ x dX- (83) 

To get the contribution to A2 at rip from the particles 
in the shell, we multiply (|83|l by A2[a,x,iVp ~ Vb)], with 
A2icr,x,v) given by 

Proceeding as before we get the analog of (|59p : 



2 rvp 



aTm 



3^?'(%) Jo V+pM^ 

dxF^ W,X, iVp - Vs)] 



(84) 



(85) 



The collision time, r, will be given by 

(1+p) 

87r In A(0.24)(?2yl^/'' (ry,) (2 + pY^/^ ' 

We can now get an expression for p by using (|85p in 
84ll and dividing by T^{rjp): 



PiVp 



Vp 



256nHnAAli\r!p) Jo 



dj7s X 



15/2 



1'2+pM] / JO 

with F given by (f77)l 



dxF'[a,x,{Vp'V^)], (86) 



15.1 Calculation of p and Am 

As before, we set InA = 10 for definiteness. Simultaneous 
numerical solution of (|8ip and (|86p is straightforward, but 
we need a starting value for Am- We choose Am = 2.5 x 
10^* eV* ■ m*. This is approximately equal to the observed 
value today, p becomes equal to unity around rj^^ — 480. 

This value of r^cw is significantly larger than the previous 
value of 11 when cr = 0, but is still within a few orders 
of magnitude of unity. The most important feature of this 
calculation is that it shows p{ri) does continue to rise, even 
when cr 7^ 0. 



16 APPLICATION TO THE MANNHEIM 
MODEL 



In a series of papers (see iMannheimI (|2006| ) and references 
given there), Mannheim has developed a cosmology based 
on conformal gravity. The underlying geometry is a FRW 
space with negative curvature. The expansion parameter, 
Riri), has a minimum, non-zero value; we will assume this 
represents the initial size of the bubble. 

A distinguishing aspect of this theory is that it contains 
no intrinsic mass scale, unlike conventional cosmology where 
the Planck mass provides a fundamental scale. It is therefore 
natural to ask whether, in Mannheim's model, the steady in- 
crease in the ratio p could eventually trigger a phase change 
that would generate particle masses. We are not advocating 
the Mannheim model in this paper; it must still meet some 
serious challenges before it can be regarded as a serious ri- 
val to the conventional model. We are simply using it as 
an example of the way the non-thermal component of the 
effective mass can play a crucial role. 

The mechanism we have in m ind is the Coleman- 
Weinberg (CW) transition (jColeman fc Weinberd 
1973ll. as discussed in the context of cosmology in 
Narlikar fc PadmanabhanI ([1986 ') . Figure [10] shows the 



with Fo given by f77\ . 



effective potential, Vch, for a typical field theory at fi- 
nite temperature. (The formulae used are taken from 
iNarlikar fc PadmanabhanI l|l986l )'). Initially, Vctt is de- 
scribed by the upper curve, but as the effective mass of the 
field quanta increases, the second minimum will move lower, 
and eventually, when p ~ 1, it will cross the horizontal 
axis. At this point a CW transition can take place and 
normal masses will appear. These will have magnitude 
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Figure 10. Effective potential, Vcfl, as a function of <I>. As 
mnt/"ith increases, the form of Vcff changes from the upper curve 
to the lower, at which point a Coleman- Weinberg transition can 
occur. 

By adjusting the parameters of the model we can ar- 
range for sufficient conformal time to elapse from the forma- 
tion of the bubble to the CW transition, assumed to occur 
at the weak interaction scale. The initial temperature of the 
bubble, Tmax, is then found to be very high, and we have to 
look for an explanation of why particle masses are so very 
small compared to this temperature. (This problem is analo- 
gous to the hierarchy problem in ordinary cosmology, where 
the Planck mass is so much greater than particle masses.) 

In the Mannheim model, the explanation is straightfor- 
ward. The temperature of the CW transition is found to de- 
pend exponentially on the conformal time, so a quite modest 
value of rj will produce a very small value of Tew /Tmax. 



17 CONCLUSION 

We have studied the propagation of the vector potential in 
the plasma of the early Universe, assuming a FRW space 
of negative curvature, and conclude that (A^) will slowly 
acquire a non-thermal part, (A^t), due to distant matter. 
Even though this is a direct effect of the vector potential, 
and cannot be attributed to fields, the theory is nonethe- 
less gauge invariant. This effect is th erefore in the same 
category as the Aharonov-B ohm effect l|Aharanov fc Bohml 
[l959„; Tonomura et al.lll982l ). 

Further research will be needed to determine whether 
this increase of (A^t) is sufficient to cause a CW transition 
and generate a mass scale. 

Several additional questions are raised by the present 
paper, among them the following: 

(i) Will a CW transition necessarily take place as (A^)nt 
approaches (A^)th ? 

(ii) Is the tail of A important only in the early Universe, 
or are there other places, regions of high gravitational fields, 
where its effects can be observed even now? 



(iii) How does the real complicated early plasma deter- 
mine such things as the collision time? 

(iv) Can we expect A to propagate over cosmological dis- 
tances so that the CW transition takes place within the 
available conformal time? We have taken findings from or- 
dinary plasma theory and used them in the very different 
circumstances of the early Universe. 

These considerations are beyond the scope of this pa- 
per, which is solely concerned with the interplay, in a simple 
model, of field theory (density matrix, the Lagrangian for a 
scalar field, the Coleman- Weinberg transition) and the clas- 
sical equations of propagation of the ordinary vector poten- 
tial in a FRW space of negative curvature. 

A well-known text (|Peskin &: SchroededllQQSh suggests 
a connection between large and small scales similar to the 
one explored here. After surveying the difficulties faced by 
current theories of the mass scale, the authors write: "... it 
may be that the overall scale of energy-momentum is gen- 
uinely ambiguous and is set by a cosmological boundary con- 
dition." We have presented a mechanism for such a connec- 
tion. It is based on the familiar electromagnetic interaction, 
and nothing radically new seems to be required. 
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APPENDIX A: NORMALIZATION WITH 
CONDUCTIVITY 

In this appendix we will work in the usual spherical polar 
coordinates, and make connection with Riemannian coordi- 
nates when necessary. 

Suppose we have a dipole at the origin, oscillating with 
time dependence exp(— itjt) in the z direction. The surround- 
ing medium is of uniform conductivity, ct, so that the cur- 
rent density, j, is given by j = ctE. We will analyze this 
system by imagining a small sphere of radius r\ cut out of 
the medium surrounding the dipole. Induced currents flow- 
ing in the medium will cause surface charges to appear on 
the sphere, and the total dipole moment will be the sum 
of the original dipole moment and that due to the induced 
charges. We assume the permittivity and magnetic suscepti- 
bility are essentially unity, so D = E and B = H. In such a 
system V • j = follows from Maxwell's equations, so there 
are no volume charges in the medium. 

Denote by Dtrue = -D(aj) exp(— io^f) the dipole moment 
at the center of the sphere, where D{uj) is the true dipole 
strength at angular frequency uj. The induced dipole mo- 
ment due to the surface charges is Dind, so the total dipole 
moment is Dtot = -Dtruc + -Dind- 

Just outside the sphere the electrostatic potential is 
given by 



1^ = AotPi(cos6»)/r?, 



(Al) 



where Pi (cos 

by 



■ cos 6. The radial component of E is given 
= 2 Act Pi (cos 61) /r?. (A2) 



The surface charge density, s, obeys the relation 



As 



(A3) 



where jr is evaluated just outside the sphere. This gives 



= soPi[cos6), 



where 

— 2iCT_Dtot 

So = 3 ■ 



(A4) 



(A5) 



The induced dipole moment, Dind, is then given by an 
integral over the surface of the sphere: 



giving 

Dtot 



2-K / d6lsin6'r^[soPi(cos6l)][rJPi(cos6')] 
Jo 



(^) (1) 



Otruc + find 
3cj_Dtruc 



(A6) 



(A7) 



3io + Siria 

We set — Cz,{io,a)Ni{6)h{r)eicp{-~iijjt), where 

6*5(0;, a) is a normalizing function, A^i = dPi/dS, and h{r) 
satisfies 



A^h 2 Ah 2h 2, „ 

-TT ^ \ 7+K h — . 

dr^ r ar 



(A8) 



with — uj{ui + 47ri(j). We choose the solution that repre- 
sents outgoing waves, so 

h{r) = h'^^\nr) 



(A9) 



(Kr)2 



— I exp(iKr). 



Here fti^^(/tr) is the s pherical Bessel function defined in 
lAbramowitz fc Stegur] ^970 ). 

The Maxwell equation V x H = — i(ci; -f- 47ria)E then 
gives 

— 216*5(0;, ct)Pi (cos 6)h\^'' (kt) exp(— itjt) 



r(tj 4- 4nia) 
For small r this becomes: 

— 2C5(tJ,CT)Pi(cosS) exp(— iijji) 



But also, for small r, we have HA2[) . so 
iuK^ {uj + 'iTiia)D(ijj) 



6*5(0;, ct) = 



{iuj -f 87rio") ' 
3ojK^(o; -I- ATTm)D{uj)Nx{e)h^^\nr) exp(-itjt) 



{iu} + 87ricr) 



6iojK D 

r(3t<j + 87ria) 



(AlO) 



(All) 



(A12) 



(A13) 



(A14) 



In this appendix we have used ordinary polar coordi- 
nates in fiat space. We need now to transform to coordinates 
of a FRW flat space, with metric 

As^ = R^{t)) [-Arf + Ax^ + % AQ^ % sa^ Q A^'') . (A15) 



To get from of (fUS)) to of (fHl). in the flat FRW 
metric (|A15I) . we first refer to [Weinberg l| 19721 ) . section 4.8, 
and multiply by x- We then convert from t\,or\ and r to x 
by multiplying by R^(r\). We also convert o; to n, k to fc and 
a to (J = R(ri)a: 



F12 



-Snk'jn + 47via)D{n)Ni{6)xh^^\kx) exp(-in77) 

P(ry)(3n + 87ria) ' ^ ^ 

We can obtain the corresponding function in curved 
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space by multiplying (|67|l by A''i(^) exp(— in77): 



Fl2 = 

C3{n,a)Ni{e)exp[i{kx - n-q)] 
2u 

We can now find C3(n, a) by matching (|A16[) with (jA17[) 
for small x'- 



(l - 2[ku + M^) . (A17) 



Cain, a) = 



3in(n + 47ri(T)_D(n) 



i?(77)(3n + 87ricr) 

With D{n) given by pop . this reduces to 

, , — 6gV(n + 47ri(T) sin(nr) 

C3(,n,o-j = — — . 

n[3n + S-Kia) 

In the limit a ^ 0, HA19|) tends to ([321), as it should 



(A18) 



(A19) 



APPENDIX B: THE MODIFIED LORENZ 
CONDITION IN CONFORMALLY 
INVARIANT THEORIES 

It is well known that although Maxwell's equations are con- 
formally invariant in 4-dimcnsional space, the Lorenz condi- 
tion is not ((Hovle fc Na rlikar 1974,, chapter 2). So a theory 
that aspires to be completely conformally invariant, such 
as the theory of Mannheim discussed in section 1161 must in- 
clude a modified gauge condition if it is to describe the prop- 
agation of potentials. We note that the conformal weights 
of ff'"^! ^ S'lid Ai, are —2, 4,-1 and zero, respectively. 
The equation 

^(Vff5'"^-J'**A.)=0 (Bl) 

is therefore conformally invariant, and is a candidate for a 
gauge condition. We can linearize it by setting $*"!> equal to 
its expectation value. In the conditions of the early Universe, 
this will be proportional to T^, where T is the temperature. 
Assuming adiabatic expansion, will be proportional to 
l/R^ , so our gauge condition reduces to 

^{^g^'^A./R^)=0, (B2) 
which is the modified Lorenz condition, (1441). 
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